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We derive macroscopic traffic equations from specific gas-kinetic equations, dropping some of the 
assumptions and approximations made in previous papers. The resulting partial differential equa- 
tions for the vehicle density and average velocity contain a non-local interaction term which is very 
favorable for a fast and robust numerical integration, so that several thousand freeway kilometers 
can be simulated in real-time. The model parameters can be easily calibrated by means of empirical 
data. They are directly related to the quantities characterizing individual driver-vehicle behavior, 
and their optimal values have the expected order of magnitude. Therefore, they allow to investigate 
the influences of varying street and weather conditions or freeway control measures. Simulation 
results for realistic model parameters are in good agreement with the diverse non-linear dynamical 
phenomena observed in freeway traffic. 
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I. INTRODUCTION 

Recently, traffic dynamics has become interesting to 
a rapidly growing community of physicists. This is not 
only due to its practical implications for optimizing free- 
way traffic, but even more because of the observed non- 
equilibrium phase transitions and non-linear dynam- 
ical phenomena like the formation of traffic jams [fjlfl. 
stop-and-go traffic ||, and synchronized traffic [j|||7 It 
seems that all forms of congested traffic have almost 
universal properties which are largely independent of 
the initial conditions and the spatially averaged density, 
like the characteristic outflow from traffic jams of about 
1800 ±200 vehicles per kilometer and lane or their typical 
dissolution velocity of about — 15 ± 5 kilometers per hour 
0. This universality arises from the highly correlated 
state of motion produced by traffic congestions H|| . 

Whereas classical approaches focussed on reproduc- 
ing the empirically observed flow-density relation and 
the regime of unstable traffic flow, recent publications 
pointed out that it is more important to have traffic mod- 
els which are able to describe the observed spectrum of 
non- linear phenomena and their characteristic properties 
H||7|,|l(J . We think that it would be desireable to de- 
velop models that are consistent with both aspects of 
empirical data. Such models have been proposed lately, 
including cellular automata models JlT|,p[ and "micro- 
scopic " models of driver- vehicle behavior |l^,[ll| , and the 
macroscopic model discussed in this paper. 

In order to have meaningful and measurable model pa- 
rameters, we will relate our macroscopic model of freeway 
traffic to a "microscopic" model of driver vehicle behav- 
ior via a gas-kinetic derivation (cf. Sec. |0]). Derivations 
of this kind have been already proposed in a number of 
previous publications Jbi]-|l7||, but the correct treatment 
of the most interesting regime of moderate and high den- 
sities remained a problem. In (l^-|2^], the effect of ve- 
hicular space requirements has been taken into account 
by a correlation factor reflecting the increased interaction 



rate of vehicles (which can now b e de rived from simple 
and plausible arguments, cf. Sec. II B). In deriving the 
associated macroscopic equations, different approxima- 
tions have been suggested, the most harmless of which 
was a gradient expansion j^^jj. This lead to a viscos- 
ity term and some unexpected but essential high-density 
corrections containing spatial derivatives of different or- 
ders. However, the resulting partial differential equations 
were not very suitable for numerical simulations. 

Meanwhile we managed to evaluate the Boltz mann - 
like gas- kinetic interaction term exactly (cf. Sec. [II B| ). 
Since it turned out that a dynamical variance equation 
p|.p2[ is not necessary for a description of the presently 
known properties of traffic flows, we replaced it by a 
constitutive relation which — corresponding to a quasi- 
adiabatic approximation — agrees with the equilibrium 
variance. The resulting macroscopic traffic equations are 
coupled non-linear partial differential equations which 
can be represented in form of flux equations with a non- 
local and anisotropic source term E3]. For this reason, 
we can now apply various standard methods for numeri- 
cal integration. It turns out that the non-local term has 
similar smoothing properties like a viscosity term, but it 
does not change the hyperbolic character of the partial 
differential equations to a parabolic one, and it has more 
favorable properties with respect to numerical stability. 
For this reason, our model allows a robust real-time sim- 
ulation of freeway stretches up to several thousand kilo- 
meters on a usual PC. 

Compared to previous macroscopic traffic models, the 
gas-kinetic-based traffic model (GKT model) proposed in 
the following takes into account the variance of vehicle 
velocities, which is basically proportional to the square 
of average velocity (cf. Eq. but with a density- 

dependent prefactor that determines the exact form of 
the flow-density relation in equilibrium. Moreover, the 
"optimal velocity" or "dynamical equilibrium velocity" 
V c (cf. Eq. (|3l|)), towards which the average veloc- 
ity relaxes, depends not only on the local density, but 
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also on the average velocity, and, even more important, 
on the density and average velocity at an "interaction 
point" which is advanced by about the safe distance. 
Nevertheless, the equations are structurally related to, 
for example, the Kerner-Konhauser model ||, so that we 
find many similar non-linear phenomena. This includes 
the sequence of stable, linearly unstable, and metastable 
regimes the local breakdown effect pq |, the local 

cluster effect ||, and the formation of dipole layers j24|] 
at sufficiently large densities. We also obtain that, in the 
unstable traffic regime, the resulting flow-density relation 
differs from the equilibrium one (lying below the latter, 
cf. Fig. ||). In addition, we find that the outflow from 
traffic jams is independent of the initial conditions and 
the spatially averaged density (cf. Fig. Moreover, 
the dissolution velocity of traffic jams varies only a little 
with density (cf. Fig. |?]). Finally, the gas-kinetic-based 
model is able to explain the common phenomenon of sy- 
chronized congested traffic |2(| , if the inflow at on-ramps 
is taken into account. 

It turns out that our model can be easily calibrated to 
the static and dynamic properties of traffic flow data by a 



certain systematic procedure (cf. Sec. HID). All param- 



eters have a clear interpretation, since they are related 
to quantities characterizing the driver vehicle units like 
desired velocities or vehicle lengths (cf. Table . At least 
some of the quantities like the typical desired velocity or 
the average time headway are directly measurable. More- 
over, the optimal parameters obtained from a calibration 
to empirical data have the expected order of magnitude 
(cf. Table |). Therefore, the model allows to investigate 
the effect of speed limits, of a larger percentage of trucks, 
of bad weather conditions, etc. It will usually be suffi- 
cient to change the affected parameter values accordingly, 
instead of calibrating or even modifying the whole model 
for every new situation. 



II. THE MODEL 

A. Underlying Gaskinetic Equation 

Similar to the gas-kinetic derivation of macroscopic 
equations for fluids, we start by formulating a kinetic 
equation for the locally averaged dynamics of driver- 
vehicles units, which play the role of the molecules, 
here. The kinetic equation describes the evolution of the 
coarse-grained phase-space density 



p(x,v,t) 



E 



dt' / dx' \ dv'g(t -t',x- x' , v - v') 



x 5(x' -x a (t)) 6(v'-v a (t)), 



(1) 



denoting the probability density of finding, at a given 
time t, a vehicle a at position x a with velocity v a . In 
more intuitive terms, p corresponds to the spatial den- 
sity of vehicles per lane times their velocity distribu- 
tion. Since the GKT model is an effectively one-lane 



model where lane changes and overtaking are only im- 
plicitely taken into account, there is no lane index. The 
coarse graining is performed by taking local averages over 
a weighting function g(t — t',x — x 1 ,v — v') satisfying 
J dt' J dx' J dv' g(t — t',x — x' , v — v') — 1, which is local- 
ized in a microscopically large and macroscopically small 
neighborhood around x and in suitable neighborhoods 
around t and v. The particular choice of g is not relevant 
for the form of the macroscopic equations p7|] . How- 
ever, the scales At, Ax, and Av are meaningful in that 
they enter the effective relations for the higher velocity 
moments like the variance. 

The information specific to vehicular traffic is con- 
tained in the "microscopic" dynamics of individual 
driver-vehicle units. In the GKT model, we assume the 
vehicular dynamics of the form [fl3|]2"s|l 



dt 7 



faf. 



(2) 
(3) 



The first term on the right-hand side of Eq. (||) represents 
the acceleration of the driver- vehicle unit a to the desired 
velocity v° with an adaptation time of t q . On empty 
roads, this is the only acceleration term. Notice that f° 
is an intrinsic property of the driver- vehicle unit. The 
second term on the right-hand side of Eq. (|^) represents 
the braking interaction of vehicle a due to slower vehicles 
(3 in front. It depends mainly on the subjective minimum 
safe time headway T a to the car in front that driver a 
wants to keep. The details of the braking interaction will 
be discussed below. 

In general, the parameters r a , and the parameters 
of the braking interaction like T a are different for each 
individual vehicle a. This could be respected by gener- 
alizing p to a multi-dimensional phase-space density in a 
phase space spanned by the dimensions x, v, t, v , etc. 
This density would express the probability of finding at x 
a driver- vehicle unit with velocity v, whose microscopic 
parameters are r, v , etc. Paveri-Fontana applied this 
concept to the extra variable v° alone, formulating an 
equation for p(x, v, v °, t) Jf5| , which was further investi- 
gated by Helbing JL7) and Wagner et al. @. 

In the GKT model, however, we assume that all de- 
viations of the individual driving behavior from that of 
the "average driver" eventually lead to fluctuations of the 
velocity. For a nearly empty road this is obvious. Then, 
the braking term is negligible, and the distribution of 
vehicle velocities must converge to that of their desired 
velocities. In dense traffic, there are additional sources of 
velocity fluctuations. The time when vehicle a starts to 
brake in response to a slower vehicle in front depends on 
the individual safe time headway T a . Obviously, driver 
a brakes later for lower values of T a . Thus, different 
T a lead to different velocities, even if all other parame- 
ters of the vehicles (in particular the desired velocities) 
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and all initial conditions are unchanged. In addition, im- 
perfect driving behavior such as delayed acceleration or 
overbraking (i.e., braking more than necessary in a given 
situation) contributes to the velocity fluctuations Jl8j . 

If we are not interested in microscopic details, such 
fluctuations can be described in a global way by a fluctu- 
ating force in the acceleration equation (0) , in analogy to 
hydrodynamic fluctuations p9[ |. Thus, we approximate 
Eq. (|) by 



dv a 
~dt 



V 



^2 fa/3 



(4) 



where V a = (v°) := jjY, a v a and l / T = QI T a) are tne 
averaged microscopic parameters of the accelaration term 
(TV denotes the total number of vehicles). The braking 
term f a g is formulated with averaged parameters like 
T = (T a ) as well. The fluctuating forces £ a (t) obey 

(£«(*)) = 0, = 2D5 ap 8{t - t'). (5) 

The fluctuation strength D will be determined empiri- 
cally by comparing the resulting velocity variance with 
that obtained from single-vehicle data, see below. Fol- 
low traffic densities, where the interactions can be ne- 
glected, Eq. (g) with (|) is an ordinary stochastic differ- 
ential equation for v a . In the stationary limit, it leads to 
the distribution function w(v a ) = (2ttDt)~ 1 exp[— (v a — 
Vo) 2 /(2Dt)]. This means that, in the low-density limit 
of negligible interactions, the fluctuating strength D of 
the fluctuating forces (0) is related to the velocity vari- 
ance 8 = ((v a — Vp) 2 ) by 6 = Dt (fluctuation-dissipation 
relation). In Sec. II B, it will be shown that this relation 



holds for stationary traffic at all densities. 

The equations (0) with (0) represent a microscopic 
traffic model on its own. It remains, however, to specify 
the "microscopic" braking interactions f a g . In real traf- 
fic, these interactions depend in a complicated manner 
on the own velocity, and on the distances and velocities 
of the vehicles in front. To formulate reasonable assump- 
tions for the GKT model, we will use averaged quantities 
like the density p{x,t) = j ' dvp{x,v,t), which is justified 
since we want to derive a macroscopic model. As a conse- 
quence, however, a simple purely microscopic equivalent 
of the GKT model exists only in some special cases. 

Specifically, we make the following simplifying assump- 
tions: i) A driver at position x reacts to the traffic situ- 
ation at the advanced "interaction point" 



x & = x + 7(Z + Tv), 



(6) 



where I = l/p max (with maximum density /? max ) is the 
average vehicle length plus the bumper-to-bumper dis- 
tance kept in standing traffic, and 7 is an average antic- 
ipation factor with typical values between 1.0 and 3.0. 
In the limit of congested traffic, the interaction point 
is 7 vehicle positions in front of the actual vehicle po- 
sition x. Notice that, in the limit of congested traf- 
fic and for 7=1, assumption i) corresponds to clas- 
sical car- following models @(H|2*|l|]. ii) In the light 



of a mean-field ansatz, the traffic situation at the in- 
teraction point can be described by the density and the 
velocity distributions at this place, i.e., by the phase- 
space density p(x a ,v,t). iii) There is a certain percent- 
age p(p a ) = l/x(Pa) of interaction-free space that al- 
lows drivers to approach the respective car in front be- 
fore they brake. This percentage is a decreasing function 
of the average density p a = J ' dvp{x a ,v,t) at the inter- 
action point with p(0) = 1, and p(p m ax) = 0. It is de- 
termined by the condition that, in homogeneous dense 
traffic, the vehicles follow each other with a time head- 
way of T . Furthermore, we assume that the probability 
of undelayed overtaking by lane changing is given by p 
as well pij| . Notice, that the factor x — \jp can be in- 
terpreted as the increase of the interaction rate due to 
the finite space requirements and positional correlations 
of vehicles, compared to point-like objects Jl8|- p0| , p8[ . iv) 
If a driver is faster than the velocity v a at the interac- 
tion point, and if he cannot overtake by lane changing 
[what happens with probability (1 — p)], he reduces his 
velocity abruptly to v a as soon as the distance to the in- 
teraction point (moving with velocity v a ) has decreased 
by his share A:Ef ree = p/ Pa of the interaction-free space 
p2j. (1/ p & is the average center-to-center distance be- 
tween two vehicles at x a ). 

Now, we formulate the kinetic equation underlying the 
GKT model. Taking the time derivative of the definition 
(0) of the phase-space density, and inserting the micro- 
scopic equations (0) and (0), gives, by partial integration, 
the kinetic evolution equation for the phase-space density 



dp d . d 



-Vb 

p — 







dv 



(7) 



where the interaction term has the general form 



/ lnt = p^^2^2 dt' dx' dv'g(t-t',x-x',v-v') 

xj aa 5[x'-x a (t)}5[v'-v a (t)}. (8) 

The four assumptions for the microscopic braking in- 
teractions f a a directly result in a Boltzmann-likc inter- 
action with a density-dependent prefactor P(p): 



^(pf i nt)=P(p)I(x,V,t) 



(9) 



with 



l(x, v,t) = j dv'(v - v)p(x, v', t)p(x a , v, t) 

v'>v 

dv'{v - v')p{x, v, t)p(x a , v', t). (10) 

The first term of ( |l0| ) describes the increase of the phase- 
space density p(x, v, t) due to the deceleration of faster 
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vehicles with velocity v' > v which cannot overtake ve- 
hicles at x a driving with velocity v, whereas the second 
term delineates the decrease of the phase-space density 
due to decelerations of vehicles driving with v which can- 
not overtake slower vehicles at x a driving with v' < v. 
The prefactor 



P(p) = (l-p) x =l-l 
P 



(11) 



is proportional to the probability (1—p) that one cannot 
immediately overtake a slower vehicle, and to the corre- 
lation factor x = 1/p describing the increased interaction 
rate due to vehicular space requirements. 

In summary, the kinetic phase-space equation upon 
which the GKT model is based, is given by 



dp d . d 



-Vo-v 
P 



= (i-i 

,P 



dv'(v' - v)p(x, v', t)p(x a , v, t) 



v'~>v 



dv'(v - v')p(x, v, t)p(x a , v', t) 



V 1 <v 



dv 2 



(pD). 



(12) 



For 7 = 1 and for dense traffic, the underlying mi- 
croscopic dynamics is that of a microscopic, stochas- 
tic car-following model. In this case, traffic behaves 
like a one-dimensional gas of inelastic hard "vehicular 
molecules" with anisotropic interactions whose effective 
sizes vary with the local density such that there is a space 
Aa:f re e(p a ) between the molecules. 



B. Derivation of the Basic Equations 

Following the standard procedure summarized in 
Refs. [^,|l8), we derive from the kinetic equation (O) 
macroscopic equations for the lowest velocity moments. 
In particular, we are interested in the dynamics of the 
macroscopic vehicle density p(x, t) per lane and the av- 
erage velocity V(x, t) defined by 



p(x,t) — / dv p(x,v,t), 



(13) 



V{x,t) = (v) = p- 1 j dv vp(x,v,t). (14) 
o 

As usual, one obtains an infinite hierarchy of equations 
where that for the nth moment depends on the (n + l)st 
moment. In particular, the macroscopic density equation 
depends on V, and the macroscopic equation for V on the 
variance 



6(x, t) = {{v ~ V) ) = p- 1 J dv (v- Vyp{x, v, t). (15) 

In the GKT model, we close the hierarchy by two assump- 
tions. First, we assume that the variance 9 is a function 
of density and average velocity. Second, we asume that 
the phase-space density is locally associated with a Gaus- 
sian velocity distribution 

e -[v-V(x,t)] 2 /[2e(x,t)] 

p(x,v,t) = p(x,t) . (16) 

^ZTTO^X, t) 

This is well compatible with empirical velocity distribu- 
tions obtained from single- vehicle data |2^j3^| , at least, 
if the percentage of trucks is negligible p4fl . Ansatz ( |l6| ) 
is also consistent with the fluctuating force @ in the mi- 
croscopic Eq. (^). A more general ansatz taking into 
account small deviations from local equilibrium can be 
found in E3l. 



Multiplying the phase-space equation (12) with 1 or v, 
respectively, and integrating over v leads, after straight- 
forward but lengthy calculations, to 



dp , d( P V) 



d_ 

di 



V 



dx 



at 



v = 



dx 

ld{pff) 
P 



0, 



(17) 



Vn-V 



dx 



-B(S V ), (18) 



where we used the notation f a (x,t) = f(x a ,t) with 
/ € {p,V,9}. It turned out, that the approximation of 
the sum (9 + 9 a )/2 » 9 leads only to negligible quanti- 
tative changes [ p5[ . On the other hand, the approxima- 
tion simplifies the velocity equation considerably, so we 
will adopt it henceforth. The monotonically increasing 
macroscopic interaction term 



B(5 V 



e-4/2 



/2tt 



dy ■ 



-y 2 /2 



(19) 



describes the dependence of the braking interaction 
on the dimensionless velocity difference Sy — (V — 
V a ) I \/9 + 9 a . For 7 = 1, the macroscopic interaction 
term can be easily understood by the underlying micro- 
scopic dynamics of the GKT model. If a vehicle at loca- 
tion x with velocity v is faster than one at x a with velocity 
v a (i.e. Sv = v — v a > 0), it approaches the car in front 
within the time At = Axf TCC /Sv, where Axf ree = Vl Pa 
is the average interaction-free distance Axf rec = p/ p a of 
a car. Then, if it cannot overtake immediately, which 
would happen with probability (1—p), it abruptly reduces 
the velocity by Sv. The resulting ensemble-averaged de- 
celeration is 



(Sv/At) =-V—E). pa 



d(6v) (6v) 2 w(6v) . (20) 
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If v and Da. are uncorrelated and Gaussian distributed [cf. 
Eq. (|l6|)], with expectation values V, V a and variances 

0, (9 a , respectively, the distribution function w(Sv) of the 
velocity difference is also a Gaussian, with expectation 
value V—V a = ^J6 + 9 a 8y and variance (9+9 a ). Evaluat- 
ing integral © yields (Sv/At) = -\Pp a {9 + 9 a )B(S v ), 

1. e., the macroscopic braking term in Eq. (|l8[). 

The assumption of a Gaussian velocity distribution 
alone would close the system after the variance equation, 
which can also be derived from the kinetic equation ( |l2] ) 
p2| , p3| . Since a dynamic variance equation is not nec- 
essary for the description of known traffic instabilities, 
we close the system already after the velocity equation 
and assume for the variance the local equilibrium value 
8 = Dt of the variance equation 22 2^) . Notice that this 
relation for the variance is the same relation as derived 
in Sec. [II A| for low densities. 

To complete the derivation of the GKT equations, we 
have to specify the "constitutive relation" for the quasi- 
adiabatically eliminated variance as a function of p and 
V, and the relation for the dimensionless correlation pref- 
actor P(p). The empirical data suggest that the variance 
(and thus D) is a density-dependent fraction A(p) of the 
squared velocity, 



A(p)V 2 



(21) 



and that the variance prefactor A is higher in congested 
traffic than in free traffic. For qualitative considerations, 
A can be chosen to be constant. In the following, how- 
ever, we approximate the empirical data by the Fermi 
function 



A(p) = A + AA 



tanh 



P- Pc 

Ap 



1 



(22) 



where Aq and Aq+2AA are about the variance prefactors 
for free and congested traffic, respectively, p c is of the 
order of the critical density for the transition from free 
to congested traffic, and Ap denotes the width of the 
transition. 

Now, we determine the correlation function P(p) by 
imposing that the time headways in dense, homogeneous 
traffic are given by T. Solving Eq. ( |l8| ) for stationary and 
homogeneous traffic of density p leads to the equilibrium 
velocity-density relation 



Ve(p) 




-1 ± 




with 



V 



Vo 



rpA(p)P(p) ■ 



(23) 



(24) 



This also determines the equilibrium traffic flow per lane 
by 



Q e (p) = pV e (p). 



(25) 



In the limit of high-densities, (1 — p/pmax) <C 1 (or V c <C 
Vo), this reduces to V c = V. On the other hand, time 
headways of T and average gaps of s = (X/p — 1/pmax) 
between the vehicles correspond to a velocity Vr(p) = 
s/T = {l/p - l//w)/T. Demanding V e = V T for high 
densities leads to 



P(P) 



V oP T 2 



TA(p max )(l - p//9 max ) 2 



(26) 



This expression is consistent also in the other limit of ho- 
mogeneous traffic with very low density. With Eq. (|2(j), 
the macroscopic braking term —Pp9 of Eq. ( fig} ) for ho- 
mogeneous traffic is proportional to p 2 , in accordance 
with intuition: The rate of encountering a slower vehicle 
is proportional to p. Furthermore, the probability that 
one cannot overtake immediately by lane-changing when 
a slower vehicle is encountered, is proportional to p as 
well, resulting in a proportionality to p 2 at low densities. 
The interpretation of p = (P + l)^ 1 as the percentage 
p Axf rcc of free space is also consistent for both limit- 
ing cases. For p — ► Pmaxj one has Aa;f ree -> 0. For 
p/pmax < 1, one obtains Aa; frco - Ax = V T 2 /(tA ). 
This means, vehicles on a nearly empty road would react 
to other vehicles in front (mostly by lane-changing), if 
these vehicles were closer than VqT 2 /(tAq). This is the 
net safety distance VqT times a factor TJ (Aqt) that is of 
order unity (see Table |). 



C. Discussion of the Model 

For convenience, let us summarize the equations of the 
GKT model. The traffic density and average velocity 
evolve, in absence of on- and off-ramps according to 



dp d(pV) 



dt 



dx 



= 0, 



(27) 



dt dx 



v _ l djpAV 2 ) 
p dx 

V A(p) ( p a TV 



Vn-V 



rA{ Pn 



1 - Pa/pn 



B(S V ), (28) 



where B(Sy) is given by Eq. (]I^). (For a generalization 
to cases with on- and off- ramps see Ref. p6| .) A(p) is 
the measured or assumed variance in units of the squared 
velocity, for which we use relation throughout this 
paper. 

The density equation ( p7j ) is just a one-dimensional 
continuity equation reflecting the conservation of the 
number of vehicles. Thus, the temporal change dp/dt 
of the vehicle density is just given by the negative gradi- 
ent —dQ/dx of the lane-averaged traffic flow Q — pV. 

The first term on the rhs of Eq. (E8J) is the gradient 
of the "traffic pressure" pO — pAV.lt describes the 
kinematic dispersion of the macroscopic velocity in inho- 
mogeneous traffic as a consequence of the finite velocity 
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variance. For example, the macroscopic velocity in front 
of a small vehicle cluster will increase even if no individ- 
ual vehicle accelerates, because the faster cars will leave 
the cluster behind. The kinematic dispersion also leads 
to a smooth density profile at the dissolution front be- 
tween congested traffic and an empty road, as it occurs 
when a road blockage at xq is removed at a time to. In 
this case, the first vehicles can all accelerate to their re- 
spective desired velocities. Thus, after sufficiently long 
times, the high-speed tail of the distribution of desired 
velocities translates into a distribution of vehicle posi- 
tions. 

The second term of Eq. ( p8| ) denotes the acceleration 
towards the (traffic-independent) average desired velocity 
Vq of the drivers with a relaxation time r. Individual 
variations of the desired velocity are accounted for by 
the finite velocity variance. 

The third term of Eq. ( p8|) models braking in re- 
sponse to the traffic situation at the advanced "inter- 
action point" x& = x + 7(l//9 m ax + TV). The braking 
deceleration increases coulomb-like with decreasing gap 
(l/p a — l/pmax) to the car in front (l/p a being the av- 
erage distance between successive vehicle positions and 
1/Pmax the minimum vehicle distance). In homogeneous 
dense traffic, the acceleration and braking terms com- 
pensate for each other at about the safe distance. In 
general, the deceleration tendency depends also on the 
velocity difference to the traffic at the interaction point, 
characterized by the "Boltzmann factor" B(5y)- In ho- 
mogeneous traffic, we have B(0) = 1. In the limiting case 
Sy S> (where the preceding cars are much slower), it 
follows B(8v) — 2o~y- If, in contrast, the preceding cars 
are much faster (i.e. 5y -C 0), we have B(Sy) ~ 0. That 
is, since the distance is increasing, then, the vehicle will 
not brake, even if its headway is smaller than the safe 
distance. 

The main difference with respect to other macroscopic 
traffic models is the non-local character of the braking 
term, which we obtained by derivation from realistic as- 
sumptions of driving behavior. The non-locality has very 
favorable properties with respect to the robustness of nu- 
merical integration methods and their integration speed. 
It has smoothing properties like the viscosity term used 
in the Kerner-Konhauser model ||||, but its effect is 
anisotropic. There is no smoothing in forward direction, 
which would imply that cars would react on density or 
velocity gradients of the vehicles behind them. 

The GKT model fits into the general class of macro- 
scopic traffic models ]l8|,|8| defined by the continuity 
equation ( p7| ) and the velocity equation 

d d \ 1 dV 1 - 

(29) 

In the GKT model, the "traffic pressure" V is given by 
V = p9 = P A(p)V 2 , (30) 



and the "dynamical equilibrium velocity" , towards which 
the average velocity relaxes in the actual traffic situation, 
is 



V e (p,V,p a ,V a ) 



V 



A{ P ) 



2A(p max ) \1- p a /p 



P.TV 



B(S X 



(31) 



In contrast to other macroscopic models |3j^jl6|,|3j belong- 
ing to the class defined by Eq. ( f29| ) , the "dynamical equi- 
librium velocity" depends on p and V at two different 
locations, thus introducing the non-locality. 

The five parameters of the GKT model, listed in Ta- 
ble ||, are all intuitive. Three of them, Vo, T, and p ma x, 
can be directly determined by fitting the equilibrium flow 
density relation ( p5f ) of the GKT model to measured 
flow-density data, cf. Fig. |J. The desired velocity Vq 
is determined by fitting the data at low densities by a 
straight line Q(p) — Vop. The safe time headway T 
and the maximum density /O max are determined by fit- 
ting the data at high densities by a straight line cross- 
ing the abszissa at pmax and identifying this line with 
Q{p) = (1 — p/pma,x)/T. The ensuing average distance 
1/Pmax of two cars in standing traffic must be consistent 
with the average length of the vehicles plus a minimal 
bumper-to bumper distance kept, which is about 1.5 m. 
As real traffic is stable at very low and very high den- 
sities, the above procedure of comparing the measured 
data in these density ranges with the equilibrium curve 
of the GKT model is justified. At intermediate densi- 
ties, the equilibrium curve of the model li es s omewhat 

0)- 



III it will be 



above the data (Fig. |l|). However, in Sec. 
shown that homogeneous traffic is unstable in this den- 
sity range, and that the averaged dynamic traffic flow in 
the GKT model is below the equilibrium curve as well. 

The remaining parameters r and 7 can be systemati- 
cally calibrated by means of the dynamic properties. This 
will be discussed in Sec. IV. Table [j] lists the numerical 
values resulting from a fit to traffic data of the Dutch 
motorway A9. If not stated otherwise, we used these val- 
ues in the numerical simulations of Sec. Ill, referring to 
them as "standard parameter set" . 

Notice that all parameters have realistic values. In 
particular, this holds for r which, for Vq = 158 km/h, 
would have the meaning of the acceleration time from 
to 100 km/h. Furthermore, Vo/t is limited to the average 
maximum acceleration of vehicles on a free road starting 
with zero velocity. For these reasons, a relaxation time 
t ~ 35 s is sensible for freeway traffic. (For city traffic, r 
is shorter). 

The value T — 1.8 s for the safe time headway is con- 
sistent with the rule "distance (in m) should not be less 
than half the velocity (in km/h)" suggested by German 
authorities. For other data, however, we often find that 
a somewhat smaller time headway gives a better fit. 

Since the model parameters are meaningful, it is sim- 
ple to model changes of the traffic dynamics caused by 
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external effects like environmental influences. For exam- 
ple, a speed limit would be conidered by decreasing V - 
Bad weather conditions leading to more defensive driving 
would be characterized by an increased time headway T 
and a lower value of Vq (plus a reduction of 7, if there 
is heavy fog). In rush-hour traffic, it is plausible to as- 
sume a higher percentage of experienced drivers than in 
holiday traffic, which would correspond to a higher 7. 
Effects like a varying distribution of vehicle types can be 
modelled as well. For example, a higher proportion of 
trucks would lead to a decrease of Vo and p max , but also 
to an increased value of r. 

Finally, we compare the macroscopic GKT model with 
direct simulations of microscopic models of the form (|J). 
While the microscopic model is stochastic, the determin- 
istic GKT model includes the stochasticity of real traf- 
fic by the constitutive relation ( pl| ) for the velocity vari- 
ance. Therefore, the GKT describes macroscopic effects 
of fluctuations like kinematic dispersion. The additional 
information of individual fluctuations contained in mi- 
croscopic models seems not to be of practical relevance, 
since empirical traffic data are typically available as one- 
minute averages, i.e. in terms of macroscopic quantities. 

In contrast to microscopic models, the GKT model is 
an effectively one-lane model and treats overtaking and 
the associated lane-changing manoeuvres in an overall 
way. A microscopic model would need additional as- 
sumptions and new parameters for the lane-changing de- 
cisions as well as additional assumptions about the pop- 
ulation of vehicles, e.g., the distribution of desired veloci- 
ties. Moreover, the macroscopic model can be generalized 
to simulate on-ramps, off-ramps, and lane closings, sim- 
ply by adding source and sink terms to the macroscopic 
density and velocity equations ( p7| ) and ( p8| ) |25|,^6|]. In 
microscopic models this would require the treatment of 
lane changes from dead-end lanes, which is a particularly 
difficult problem. 

Finally, apart from very low densities, the numeric per- 
formance of simulations with the GKT model is far su- 
perior to corresponding microscopic simulations. This 
is achieved partly by using lookup tables for the func- 
tions A(p) and B(8v) and by applying explicit integra- 
tion schemes p7l]. In addition, the GKT model has only 
one density and velocity field variable, and its computa- 
tional speed (measured in terms of the length of the road 
sections that can be simulated real-time) is independent 
of the density and the number of lanes. 



D. Dimensionless Form of the GKT Equations 

By reformulating the GKT model in dimensionless 
variables, the number of model parameters can be re- 
duced by two. We measure times in units of r and dis- 
tances in units of l/p max by introducing t' = t/r and 
x' = Pm&xX- The dependent variables p' = p/p max , 
V = Vrp max , 6' — ^T 2 p^ ax , etc. are scaled accordingly. 



This leads to the scaled GKT equations 
dp' d(p'V') 



Of 



dx' 







(32) 



and 



p' dx' 

(o' V') 2 

P'A'(p') \ P ^ J B(Sv), (33) 

\^ Pa) 



where 



u, f \ A{p max p') 



A'(p') 



^-(Pmax) 



(34) 



is of order unity, and the Boltzmann term B(6v) (i.e., 
Eq. (jl^) in scaled variables), remains unchanged. The 
remaining dimensionless parameters are the scaled de- 
sired velocity 



V = PmnxTVo 

and the scaled cross-section 



(35) 



(36) 



in addition to the anticipation factor 7 from the unsealed 
equations. 

The parameter V^', with a numerical value of 171.4 
for the standard parameter set, has some analogies to 
the Reynolds number in the Navier-Stokes equations for 
normal fluids. Assuming that typical velocities are pro- 
portional to the desired velocities, typical densities pro- 
portional to Pmax, arid typical length scales proportional 
to l//9 m ax, this can be seen by observing that the magni- 
tude of the destabilizing advection, pressure, and braking 
terms in the unsealed velocity equation (^8|) is propor- 
tional to /O max Vg 2 , while the stabilizing relaxation term is 
proportional to Vq/t. So, the ratio between the destabi- 
lizing "kinetic" terms and the stabilizing relaxation term 
is p ropor tional to p ma xVbr = Vq. As will be shown in 
Sec. [ill B| , homogeneous traffic can become unstable, if a 
certain "critical" value V<! is exceeded. 

The scaled cross-section J", with a numerical value 
of 0.453 for the standard parameter set, gives the ratio 
between the interaction term and the "kinetic" advection 
and pressure terms. In analogy to the Prandtl number of 
thermal convection in a simple fluid heated from below 
|p8| ) it depends on the ratio of the two relevant time 
scales T and r of the system. 

As in the unsealed equations, the parameter 7 charac- 
terizes the sensitivity in the braking interactions to spa- 
tial changes of the velocity, compared to the sensitivity 
for changes of the density. 
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III. RESULTS 
A. Homogeneous Traffic 

In homogeneous and stationary traffic, the GKT equa- 
tions (|27]) and ( p8| ) reduce to p — const, and relation 
( p3| ) for the equilibrium velocity V e (p). Notice that V c (p) 
and thus the equilibrium flow Q c = pVc(p) is a function 
of the model parameters Vq, T, and p max , but does not 
depend on r and 7. 

We determined the constants in the relation ( p2] ) for 
the variance prefactor A(p) and the model parameters 
Vb, T, and p max , by fitting them to empirical data of the 
Dutch two-lane motorway A9 from Haarlem to Amster- 
dam. The empirical data are based on one-minute values 
for the number rit of passing vehicles, their average veloc- 
ity Vf, and velocity variance Q t , which were determined 
from single-vehicle data. The corresponding flow is then 
given by Qt = nt/2 per minute and lane, and the den- 
sity per lane by Qt/V t . Such sets of one-minute values 
were sampled over two periods from Monday to Friday 
(October 10, 1994 through October 14, 1994 and Octo- 
ber 31, 1994 through November 4, 1994) for the right 
and left lane at nine subsequent measuring cross-sections 
(distributed over a stretch of 8.6 km length; see Ref. |||] 
for an illustration). To obtain empirical quantities as a 
function of density, we averaged over all sets with a den- 
sity between p— Ap and p + Ap with Ap —1 vehicle/km. 

Figure |l](a) shows the variance prefactor Eq. (^2[) for 
the fitted values 

p c = 0.27p max , 
Ap = 0.05p max , 

= 0.008, and 
AA = 2.5A . (37) 

Throughout the paper, we will use these values. 

Figure |l|(b) depicts the flow-density relation for the 
fitted values of Vb, T, and p max given in Table 61 These 
values r esult ed from the systematic procedure described 
in Sec. |III D . Both the empirical variance-density rela- 
tion and the flow-density relation are well reproduced by 
introducing one fit function A(p) only. With a constant 
value for A, the sharp increase of the variance prefactor 
at a density of about 40 vehicles/km and the sharp de- 
crease of the velocity related to it could not be obtained 
by variation of Vb, T, and p max . Rather, this correlation 
is an intrinsic property of the model and follows from the 
proportionality of the braking interaction to the variance. 
Notice, that the deviation of the assumed variance pref- 
actor from the data at very low densities could be easily 
removed by a more complicated function A(p). How- 
ever, this correction would not impair the flow-density 
relation, because the interaction is negligible for these 
densities. For the same reason, the correction would not 
change the dynamics, justifying the choice of the simple 
functional dependence (E3). 



In the plots of Fig. |^, the resulting equilibrium flow 
Q c = pV e {p) is plotted as a function of the density for 
various values of the model parameters Vb, T, and p max - 
The solid lines in the plots |2|(a) - ||(d) show the result for 
the standard parameter set of Table Q. 

Figures ||(a) and (b) illustrate that, according to our 
model, traffic flow will increase with growing desired ve- 
locity Vb and with decreasing safe time headway T, as 
expected. The desired velocity Vb mainly influences the 
low-density regime, while T is relevant for high densities. 
This is also underlined by Figs. ||(a) and ||(b), which 
show the equilibrium velocity for the same parameters. 
As expected, a reduction of Vb (e.g. caused by a speed 
limit or an uphill gradient) influences traffic only at low 
densities. Since a reduction of Vb increases the stability 
of homogeneous traffic (see below), this has practical im- 
plications. A variable speed limit, which is active only 
above a certain density threshold, would increase the sta- 
bility without reducing the capacity of the road. 

The changes of traffic flow resulting from variations of 
p max are plausible as well [Fig. ^(c)]. Decreasing p max or, 
equivalently, increasing the average vehicle length (cor- 
responding to a higher percentage of trucks) reduces the 
maximum flow (capacity) of the road. Of course, a higher 
proportion of trucks also leads to a lower desired velocity. 
These combined effects are shown in Fig. [|(d). 



B. Stability of Homogeneous Traffic With Respect 
to a Localized Perturbation 

The homogeneous and stationary equilibrium solution 
investigated in Sec. Ill A is not always stable. Here, 



we consider its stability with respect to a localized per- 
turbation in the initial conditions. Specifically, we as- 
sume a dipole-like initial variation of the average density 
p according to 



p(x,0) 



cosh 



Ap 

w + 

— — cosh 
w 



x 



x 



X — Xq 



Ax 



w 



(38) 



as suggested in Ref. . The positive and negative peaks 
are positioned at xq and xq + Axq, respectively, with Axq 
= 1006.25 m, and they have the widths w + = 201.25 m, 
and w~ — 805 m, respectively. The amplitude Ap of the 
perturbation will be varied in the simulations. The initial 
flow Q(x, 0) = Q c (p(x, 0)) is assumed to be in local equi- 
librium everywhere. Our simulations showed that the 
specific shape of the perturbation is not relevant. More- 
over, localized perturbations of a different form (e.g., a 
constant density plus a Gaussian perturbation of the av- 
erage velocity) were, after a short time, transformed to 
dipole-like perturbations similar to that of Eq. fl38|). The 
simulations were carried out with explicit finite-difference 
methods El. 
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Figure |j shows the spatio-temporal evolution of the ini- 
tial perturbation ( |38| ) with Ap — 10 vehicles/km for var- 
ious initial densities p. In Fig. f|(a), it is shown that the 
perturbation dissipates, if the traffic density is sufficiently 
low. When increasing the initial density, perturbations 
eventually lead to instabilities. Depending on the den- 
sity, a single density cluster [cf. Fig. |5|(b)], or a cascade 
of traffic jams (i.e., stop-and-go traffic) [Fig. ^(c)] is trig- 
gered. If one further increases the density, we observed 
dipole-like structures similarly to those described in Ref. 
p4j . Finally, for densities above 50-55 vehicles/km, one 
reaches again a stable regime [Fig. f|(d)]. 

The underlying instability mechanism is intuitive [^J. 
In the perturbation (|38|), the positive density peak (with 
a lower velocity than in the homogeneous regions) is 
behind the negative peak (with higher velocity). This 
means, drivers in the region upstream of the perturbation 
will approach the positive peak. If the traffic density is 
sufficiently low, the vehicles can overtake by lane chang- 
ing without braking, as soon as they meet the tail of the 
perturbation. In addition, the perturbation dissolves by 
the above described effect of kinematic dispersion. As a 
result, homogeneous traffic is stable at low densities. For 
higher densities, however, a higher percentage of drivers 
approaching the density peak must brake, thereby lo- 
cally increasing the density. The increased density, in 
turn, gives a positive feedback for further braking reac- 
tions. This feedback cycle continues until the resulting 
velocity is so low that the acceleration term compensates 
the braking effects. This defines the jammed density 
P = Pjam with nonzero flow Q — Q c (Pjam)- Furthermore, 
it makes plausible that homogeneous traffic of density 
p > Pjam is stable again. 

An important criterium for realistic traffic models is, 
that the transition to these inhomogeneous states should 
be hysteretic, corresponding to a first-order phase tran- 
sition j5|. This implies that, at least in some parameter 
range, the response of the system to localized perturba- 
tion depends on the perturbation amplitude (bistability). 
In terms of real traffic, there are situations where traffic 
flow is metastable with respect to small perturbations, 
but it breaks down if the perturbations are sufficiently 
large. Figure |5| shows that we found two ranges of densi- 
ties, p 6 [p c i, Pc2] and p € [pc3, p C 4], where the transition 
is bistable. For large perturbation amplitudes, the sys- 
tem develops to a localized-cluster state if p € [p c ,i, p C 2], 
or to a dipole-like state if p € [p C 3,Pc4], while it relaxes 
back to the equilibrium state for smaller perturbations. 
Between these ranges, there is a region p 6 [p C 2,Pc3], 
where homogeneous traffic is linearly unstable, giving rise 
to cascades of density clusters ("stop-and-go waves"). 
For p < p c i and p > p C 4, traffic is stable for arbitrary 
perturbations. 

In Fig. |5|(a), we plot the minimum and maximum den- 
sities pmin and pjam that resulted from simulations of 
traffic on a circular road after a dynamic steady state 
has been reached. A large difference (pjam — Pmin) cor- 
responds to localized-cluster or stop-and-go states, while 



(Pjam — Pmin) ~ is the signature of the homogeneous 
state. At a perturbation amplitude of one vehicle per 
km (dashed lines), the perturbation can be considered 
as being linear, while at Ap = 20 vehicles/km (dotted 
lines), the perturbation amplitude reaches the order of 
the homogeneous densities defining the maximum per- 
turbation. So, we can determine from the plot the four 
critical densities. Notice that the densities pj am and p m i n 
of the developed non-linear state depend neither on the 
initial density nor on the amplitude of the perturbation. 
Figure ^ shows the basins of attraction of the two traf- 
fic states in the phase space spanned by p and Ap. For 
values of (p, Ap) corresponding to points inside the two 
curves, the final steady-state consists of density clusters 
or stop-and-go waves. Otherwise, the final state is that 
of homogeneous equilibrium. 

Now, we will show that the stability of the model is 
determined mainly by the relaxation time r and by the 
anticipation factor 7. Since the flow-density relation does 
not depend on these parameters, this means that one can 
calibrate the stability and the flow rates independently. 
Figures ||(b) and ^| show that an increased value of r 
leads both to an increased range of instability, and to in- 
creased amplitudes (pj am — Pmin) of the non-linear state. 
If t exceeds values of about 60 s, the density exceeds 
the maximum density p max in the course of the simula- 
tions, which is a signature of accidents. For r < 12 s, 
the unstable and metastable regions vanish altogether, 
and the system is globally stable for all densities and all 
perturbation amplitudes. For a fixed value 7 = 1.2 but 
otherwise arbitrary model parameters, global stability is 
reached, if the dimcnsionless parameter V ' — pmaxVor of 
the scaled equations ( |32"| ) and (|3|) satisfies Vq < V c ' = 59. 
Further simulations showed that the stability of traffic 
described by the GKT model increases also with increas- 
ing 7. Both the critical relaxation time where accidents 
occur, and the critical value for global stability in- 
crease. This is plausible since 7 desribes the anticipation 
of future velocity changes. 

It turned out that the qualitative difference between 
the non-linear state of Fig. Wb), where only one clus- 
ter is seen, and that of Fig. [|(c), where a cascade of 
stop-and-go waves appeared, can be understood by the 
linearly unstable and metastable density ranges discussed 
above. In both cases, there is a region in the wake of the 
first cluster, where the density is lower, and the velocity 
is higher than in the homogeneous region downstream. 
This gives rise to a transition layer between these regions 
serving as a small perturbation || . In Fig. ||(b) , the sys- 
tem is metastable, and this perturbation is too small to 
induce a new density cluster. In Fig. ^(c), the system is 
linearly unstable and the transition layer triggers a new 
cluster. This cluster, in turn, gives rise to a new transi- 
tion layer eventually leading to a cascade of stop-and-go 
waves. 

Another remarkable feature of Fig. |^(c) is that the 
width of the first (upstream) stop-and go wave is grow- 
ing while the other waves remain narrow. The width of 
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the cluster in Fig. |J(b) increases as well, but slower. 
This can be explained by observing that the growth rate 
of the width is given by the difference of the group ve- 
locities U g! u p = (Qjam - <2in)/(Pjam - Pin) and V g = 

(Qout - Qjam)/(Pout - Pjam) of the upstream and down- 
stream fronts (Fig. |?]). For the density p = 25 vehi- 
cles/km corresponding to Fig. |J(b), the difference is 
small, while for p = 35 vehicles/km [Fig. [|(c)], it is 
large. In all subsequent clusters of Fig. £|(c), we have 
Qin = <3out and p in = p out , so they do not grow. 

It is required for any realistic traffic model, that the 
density pj am inside localized clusters and the traffic in 
the region downstream of it is independent of the inflow 
||. This implies that the group velocity v g should be 
constant. The upper curves in Fig. || for the jam density, 
the dashed line of Fig. ^ for v g , and the dashed line of 
Fig. H for Qout show, that the GKT model essentially 
fulfills this requirement. 

Another frequently discussed problem is, to which ex- 
tent equilibrium flow-density relations of traffic models 
can be calibrated by non-equilibrium empirical data. In 
Fig. ^ we simulate the measurement of traffic data by 
recording the velocity and density at several fixed loca- 
tions ("measuring cross-sections") over a certain period 
of time. To incorporate the distribution of densities en- 
countered in real traffic, we run simulations for various 
densities. In each run, the simulation time was propor- 
tional to the occurence probability of traffic densities in 
the course of the day, for which we assumed a linear de- 
crease from a maximum value at p = to zero for p 
70 vehicles/km. The dots in Figure [|show the flows and 
densities "measured" at the cross-sections for all simu- 
lation runs put together. In addition, the solid line in 
Fig. ^ shows the dynamic average of the flow for runs 
at a given density. Both figures show that, in the case 
of unstable traffic, the flow of the "measured" dynamic 
flow-density relation (broken line) is considerably lower 
than the equilibrium flow. From this it follows that, in re- 
gions of unstable traffic, one cannot calibrate equilibrium 
flows of traffic models to empirical flow-density relations. 

C. Fronts Between Congested and Free Traffic 

The realistic description of shock fronts in traffic is a 
particularly difficult problem, as pointed out in Ref. . 
Therefore, in this section we will investigate how fronts 
between two different states of traffic, especially between 
free and congested traffic, evolve in the GKT model. We 
model such fronts by initial conditions containing dis- 
continuities in the fields p and V . In particular, we con- 
sider shocks at x sn0 ck between two homogeneous regions 
with densities pi and p2 at the left (upstream) and right 
(downstream) sides. Initial conditions with pi > p2 in- 
clude situations of dissolving jams. The extreme case 
Pi = Pmax corresponds to vehicles starting from stand- 
still after a road blockage is removed. Initial conditions 



with pi < p2 include situations, where free traffic flow 
meets a queue of nearly standing vehicles. The simula- 
tions were carried out with the standard parameter set of 
Table |. At the upstream boundary (x — km), we chose 
the fixed (Dirichlet) boundary conditions p(0, t) — pi and 
Q(0,t) = Qe(pi) to model a constant inflow consistent 
with the initial conditions. At the downstream boundary 
(i ma x = 40 km) , we simulated an unperturbed outflow by 
the homogeneous von Neumann conditions d x p(x max ,t) 
= d x Q(x max ,t) — 0. Figure [l0|(a) and |l(](b) show the 
spatio-temporal development of the density and the flow 
for pi < pi (upstream jam-fronts). 

One can see that the shape of the backwards moving 
front does not change in time [cf. Fig ji~0|(a )1 and that 
there is no region of negative velocity [Fig.|10|(b)]. This is 
achieved by the non-local interaction term in the veloc- 
ity equation, while a viscosity term (~ d 2 V/dx 2 ) would 
make the front smoother with increasing time. 

Figures |ll|(a) and |ll](b) show the spatio-temporal evo- 
lution of the density and flow for downstream fronts with 
Pi 3> p2- This corresponds to a dissolving jam (for exam- 
ple, after an accident has been cleared) with an outflow 
to a nearly empty road section. Due to the kinematic 
dispersion term (l/p)d(p9)/dx in the velocity equation, 
the forward moving front is somewhat smoothed in the 
course of time. The finite velocity variance implies that, 
after some time, the faster cars are found in a wider dis- 
tance from the jam front. 

As a remarkable fact it should be mentioned that, 
although the equilibrium point of the outflow (p ou t, 
Oe(pout)) is the result of a dynamic process, the out- 
flow Qout is nearly constant over a wide range of pi [cf. 
Fig.|lJ(a)]. This agrees with empirical observations where 
it has been found that the outflow of very different forms 
of congested traffic (including the dissolution of queued 
city traffic after a traffic light turns green) is nearly a 
"constant of traffic" 

This generalizes the above mentioned requirement of 
a constant outflow from clusters (whose jam density is 
determined by the dynamics) to forms of congestion that 
are a result of the initial and boundary conditions. 

However, the outflow Q ou t varies with the model pa- 
rameters. It increases with growing 7 [cf. Fig.|l2](b)], 
because an acceleration tendency is already recognized 
in a larger distance from the jam front. The outflow 
decreases drastically with increasing T [cf. Figjl2|(c)], 
because T determines the time headway between two fol- 
lowing vehicles or, equivalently, the inverse of the flow. 
Furthermore, the outflow is diminished with increased r 
[cf . Figp^(d)] , because increased values of r correspond 
to lower accelerations and thus to more inert vehicles. 



D. Method of Parameter Calibration 

While the model parameters Vq, T, and p max influence 
the equilibrium flow-density (or velocity-density) relation 
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(Figs. H and ||) , the parameters r (Fig. |J) and 7 influence 
only the stability behavior. This enables an effective cal- 
ibration of the GKT model to concrete traffic situations. 

First, Vo is determined as the low-density limit of 
the experimental velocity-density relation (which, in this 
limit, is also the equilibrium relation), and p max is de- 
termined by the average length of the vehicles and as- 
suming a reasonable bumper-to bumper distance of, e.g., 
1.5 m in standing traffic. Then, T is calibrated by the ob- 
served maximum flows, by the outflows from stop-and-go 
waves, or by the flow resulting from standing traffic after 
a red light turns green or an obstacle is removed |p|,p3[. 
Afterwards, one calibrates r and 7 by the observed sta- 
bility behavior (Fig. |J) and by the shape and width of 
the downstream and upstream fronts connecting free and 
congested states (Figs. |l^ and p] ). 

Since r and 7 weakly influence the flows [Figs. |l2|(d) 
and [l2|(b)], the calibration of T, r, and 7 is repeated 
recursively until convergence is obtained. 

Applying this procedure to the single-vehicle data of 
the Dutch motorway A9 leads to the parameter set shown 
in Table Q. Notice that, because of their immediate intu- 
itive meaning, the plausible range of values for the model 
parameters is rather restricted. One can argue that rea- 
sonable time headways are in the range T £ [1.0s, 2.5s], 
and that initial accelerations a max = Vq/t are in the 
range a max £ [1 m/s 2 , 4 m/s 2 ] corresponding (e.g. for V 
= 144 km/h) to r £ [10 s, 40 s] . Finally, the minimum an- 
ticipation of traffic is to the car in front, implying 7 > 1. 

IV. SUMMARY AND CONCLUSIONS 

We have proposed a macroscopic gas-kinetic-based 
traffic model (GKT model) that was derived from a mi- 
croscopic model of vehicle dynamics. The assumed fluc- 
tuations of vehicle acceleration implied a velocity distri- 
bution of finite variance which is governed by a kinetic 
equation related to that used in kinetic gas theory. In 
the resulting GKT equations, the velocity variance enters 
both the braking interactions and the smoothing effect of 
the "kinematic dispersion" . 

In contrast to gas-kinetic-based models proposed ear- 
lier jl^j2l]] , we could now derive the correlation factor re- 
flecting the increased interaction rate of vehicles at high 
densities by simple and plausible arguments. Further- 
more, we replaced the dynamic variance by a "constitu- 
tive relation" obtained from single-vehicle data, thereby 
considerably simplifying the description. The resulting 
relation of the flow as a function of the density agreed 
well with empirical data. Since the form of this function 
is determined by the constitutive relation, this supports 
the assumption of the model that braking interactions 
are mainly caused by a finite velocity variance. 

Because of its derivation from physical assumptions, 
all model parameters of the GKT model have an intu- 
itive meaning and can be either directly measured or cal- 



ibrated to real traffic data. It is straighforward to model 
the effects of, e.g., speed limits, varying road conditions 
(e.g., gradients), or different driving behavior. Moreover, 
we proposed a systematic calibration procedure. 

In contrast to other macroscopic traffic models, the 
GKT model is non-local. While the non-locality has sim- 
ilar smoothing effects as a diffusion or viscosity term, it 
leads to a more favourable numerical stability behavior. 
The model belongs to the class of effective one-lane mod- 
els, so multi-lane aspects are treated only in a global way. 
Nevertheless, it is straightforward to include ramps p6| , 
or lane blockages (e.g., due to road works or accidents). 

Investigations of instabilities arising from a localized 
perturbation of homogeneous traffic showed qualitatively 
the same scenario as has been found for the model of 
Kerner and Konhauser (|J2^]. For sufficiently high val- 
ues of r, homogeneous traffic at intermediate densities 
becomes unstable with respect to perturbations, leading 
to non-linear states like localized clusters, stop-and-go 
waves, or dipole-like structures, while homogeneous traf- 
fic flow is stable for high and low densities. The GKT 
model satisfies two requirements which should hold for 
any realistic model ||. First, the transition to localized 
clusters is hysteretic. Second, the outflow region of the 
clusters is nearly independent of the homogeneous den- 
sity. In contrast to most other models, however, the traf- 
fic flow within the clusters is finite and can be influenced 
by the model parameters r and 7. This enables the sim- 
ulation of "synchronized" traffic states |2(J that turned 
out to be the most frequent form of congested traffic [Q . 

Only for high densities and for 7=1, the interac- 
tion term of the microscopic equations (Q) underlying 
the GKT model can be written in terms of a simple 
car-following model. For this case, we performed sim- 
ulations with a smoothed version of Eq. (|j) (containing 
no abrupt velocity changes in the interaction term) and 
found a nearly quantitative agreement in describing col- 
lective states like clusters or stop-and go waves. 

The high numerical stability of the GKT model also 
allowed to treat realistic boundary conditions (instead of 
the periodic boundary conditions used in most previous 
publications), and to simulate discontinuous fronts be- 
tween homogeneous low-density and high-density states. 
Such fronts correspond to the formation or dissolution 
of jammed traffic that is caused by initial or boundary 
conditions rather than by dynamic instabilities. 

Remarkably, the outflow Q on t from jammed regions 
was nearly the same as the outflow from localized clus- 
ters, regardless of the density of the jams, including even 
standstill traffic (related, e.g., to the dissolution of a 
queue behind a traffic light turning green) . The observa- 
tion that the outflow from arbitrary kinds of congested 
traffic is a "universal" constant of traffic dynamics, and 
also its numerical value of about 1800 vehicle s per hour, 
agrees well with observations of real traffic |8p4|. Be- 
sides the constant group velocity of localized clusters, 
this universal outflow can be considered as an additional 
requirement for realistic traffic models. 
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Parameter 


Symbol 


Typical Value 


Desired Velocity 


Vo 


110 km/h 


Maximum Density 


Pmax 


160 vehicles/km 


Acceleration Relaxation Time 


T 


35 s 


Time Headway 


T 


1.8 s 


Anticipation Factor 


7 


1.2 



TABLE I. Typical parameter values of the GKT model 
used for the simulations throughout this paper. The val- 
ues were obtained by calibration of the model parameters to 
Dutch freeway data. 
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FIG. 1. Comparison of (a) the density-dependent rela- 
tive variance in units of the squared average velocity, and (b) 
the equilibrium flow-density relation (^H|) in the GKT model 
(solid lines) with empirical data (crosses) . The empirical data 
were obtained from single- vehicle data of the Dutch motorway 
A9 by averaging over one-minute intervals (see main text). 
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FIG. 2. Equilibrium flow-density relations of the GKT 
model. The diagrams (a) through (c) show the variation with 
the model parameters Vb, T, and p max . In each diagram, 
the solid lines correspond to the standard parameter set dis- 
played in Table i. Diagram (d) shows parameter combina- 
tions as one would expect for, e.g., an increased percentage of 
trucks (Vb = 80km/h, T = 1.8 s, p max = 140 vehicles/km), 
or a decreased percentage (Vo = 140km/h, T = 1.4 s, 
Pmax = 180 vehicles/km). 



FIG. 3. Variation of the equilibrium velocity-density rela- 
tion of the GKT model with (a) the desired velocity Vb, and 
(b) the safe time headway T. The other parameters have the 
standard values displayed in Table [l| 
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FIG. 4. Spatio-temporal evolution of traffic on a 
uni-directional ring of circumference 10 km, starting with ini- 
tially homogeneous traffic (of density p) to which a localized 
perturbation of amplitude Ap = 10 vehicles/km is added in 
accordance with Eq. (^). (a) is for p — 15 vehicles/km (lin- 
early stable), (b) and (c) are for unstable traffic (p = 25 vehi- 
cles/km and 35 vehicles/km), and (d) is for stable congested 
traffic (55 vehicles/km). The model parameters are given by 
the standard set displayed in Table |. 
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FIG. 5. Stability diagram for perturbations of homoge- 
neous traffic of the form (Bq) in a ring of circumference 10 km. 
Both diagrams show the maximum and minimum densities 
Pjam and pout as a function of the average density p, measured 
after a dynamical equilibrium was reached. The unstable traf- 
fic regime corresponds to the density range where the jam am- 
plitude (pjam — Pout) is large (rectangle-like shaped regions). 
Diagram (a) shows the dependence of the stability diagram 
on the perturbation amplitude Ap. One can clearly see two 
density ranges [p c i, Pc2] with p c i = 21 vehicles/km and p C 2 = 
24 vehicles/km, and [p C 3,p C 4] with p C 3 = 51 vehicles/km and 
p C 4 = 55 vehicles/km, where traffic is non- linearly stable, i.e., 
stable for small perturbations, but unstable for large pertur- 
bations. In the range [p C 2,Pc3], homogeneous traffic is unsta- 
ble for arbitrary perturbation amplitudes. Diagram (b) shows 
the stability diagram for various relaxation times r and a per- 
turbation amplitude of Ap = 1 vehicle/km. 
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FIG. 6. Critical perturbation amplitudes Ap c of localized 
perturbations in the metastable (non-linearly unstable) den- 
sity regime for various values of r. The critical amplitude 
is the minimal amplitude that can cause traffic jam forma- 
tion. Smaller amplitudes are damped out in the course of 
time. Between the two metastable density regions arbitrar- 
ily small perturbation amplitudes will cause the formation of 
traffic jams. At low and high densities, inhomogeneities of 
traffic flow tend to disappear. 
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FIG. 7. (a) Group velocity of the upstream front (solid) 
and downstream front (dashed) of the first stop-and-go wave 
displayed in Fig. The propagation of the fronts is calcu- 
lated during a period, where the density of traffic jams does 
not grow anymore (t > 4 min) , but the dynamics does not yet 
depend on the boundary conditions (t < 15 min). Notice that 
the propagation velocity at low densities is positive, but slower 
than the average vehicle velocity. In the instability region, 
the negative propagation velocity of the downstream front de- 
pends only weakly on the initial density, and its magnitude is 
well compatible with empirical data. The larger propagation 
velocity of the upstream front comes from the growing jam 
length, (b) Group velocity of the upstream and downstream 
fronts on a circular road. Here, the length of traffic jams sta- 
bilizes after some time, and their upstream fronts move with 
the same velocities as their downstream fronts. The magni- 
tude of the initial perturbation (Ap = 10 vehicles/km) has 
been chosen larger than in diagram (a) (Ap =1 vehicle/km), 
resulting in a larger region of negative group velocities. 
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FIG. 8. Characteristic flows in the fully developed 
stop-and-go traffic corresponding to Fig. |B|(b), which results 
on a circular road from locally perturbed traffic of average 
density p. Depicted are, as a function of p, the flows Qjam in 
the jammed regions, the outflows Qout from jams, and the av- 
erage flows. For comparison, the equilibrium flow Q c — pV c 
with V e from Eq. (|23|) is also shown. Notice that, in the 
unstable range, the average dynamic flow is lower than the 
equilibrium flow. 
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FIG. 9. Simulated flow-density relation as superposition 
of data obtained at five equally distributed "detectors" on 
a circular road of length 10km (dots). In order to take 
into account the variation of traffic densities in the course 
of the day, the average density p was varied in the range 
[10 vehicles/km, 70 vehicles/km] in steps of 2 veh/km. The 
total duration of simulations for the homogeneous density p 
was 70min[l — p/(70 vehicles/km)]. Note that, in the unsta- 
ble regime, the dynamical flow-density values tend to lie below 
the equilibrium relation (dashed line). 
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FIG. 10. Simulation of an upstream front with initial 
densities of pi = 15 vehicles/km and p2 = 140 vehicles/km. 
Shown is the spatio-temporal evolution of (a) the density, and 
(b) the flow. In (b), the direction of the space and time axes 
is reversed for illustrative reasons. 
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FIG. 11. Simulations of downstream- fronts with Pj am = 140 
vehicles/km. Shown is the spatio-temporal development of (a) 
the density, and (b) the flow. Free boundary conditions were 
used on both sides. 
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FIG. 12. Outflow Qout from a congested traffic state of 
density pi in dependence of (a) pi, (b) 7, (c) T, and (d) r. 
The simulated traffic situation is that of Fig. nil The outflows 
were determined after a transient time of 30 minutes. 



CD 

CD 
> 



2300 



2000 



1700 



6 1400 




80 100 120 140 
Pl (vehicles/km) 



CD 

CD 
> 



2300 



2000 



1700 



6 1400 




0.8 1 



1.2 

y 



1.4 1.6 



CD 

CD 
> 



2300 



2000 



1700 



6 1400 



g 2300 
*§ 2000 
sb 1700 



O 1400 




20 25 30 35 
x(s) 



40 



